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Abstract. The aim of this article is to prove asymptotic shape theorems 
for the contact process in stationary random environment. These theorems 
generalize known results for the classical contact process. In particular, if Ht 
denotes the set of already occupied sites at time t, we show that for almost ev- 
ery environment, when the contact process survives, the set Ht/t almost surely 
converges to a compact set that only depends on the law of the environment. 
To this aim, we prove a new almost subadditive ergodic theorem, croissance 
aleatoire, processus de contact, environnement aleatoire, theoreme ergodique 
presque sous-additif, theoreme de forme asymptotique. 
60K35,82B43. 



1. Introduction 

The aim of this paper is to obtain an asymptotic shape theorem for the contact 
process in random environment on 7L d . The environment is here given by a collec- 
tion (X e )eeE d °f random variables indexed by the set E d of edges of the lattice 1 d : 
the random variable A e gives the infection rate on edge e, while each site becomes 
healthy at rate 1. We assume that the law of (A e ) e6E d is stationary and ergodic. 

Our main result is the following: if we assume that the minimal value taken by 
the (A e ) e gE<i is above A c (Z d ) - the critical parameter for the ordinary contact process 
on % d - then there exists a norm fi on R d such that for almost every environment 
A = (X e ) eGE d, the set H t of points already infected before time t satisfies: 

Pa (3T > t > T (1 - e)tA^ c H t C (1 + e)tA„) = 1, 

where is the unit ball for /i and Pa is the law of the contact process in the 
environment A, conditioned to survive. The growth of the contact process in ran- 
dom environment is thus asymptotically linear in time, and governed by an shape 
theorem as in the case of the classical contact process on Z d . 

Until now, most of the work devoted to the study of the contact process in 
random environment focuses on determining conditions for its survival (Liggett |29j . 
Andjel [3], Newman and Volchan [31]) or its extinction (Klein [27]) . They also 
mainly deal with the case of dimension d = 1. Concerning the linear growth, 
Bramson, Durrett and Schonmann [6] show that a random environment can gives 
birth to a sublinear growth. On the contrary, they conjecture that the growth 
should be of linear order for d > 2 as soon as the survival is possible, and that an 
asymptotic shape result should hold. 

For the classical contact process, the proof of the shape result mainly falls in two 
parts: 
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• The result is first proved for large values of the infection rate A by Durrett 
and Griffeath [16] in 1982. They first obtain, for large A, estimates essen- 
tially implying that the growth is of linear order, and then they get the 
shape result with (almost) subadditive techniques. 

• Later, Bezuidenhout and Grimmett [4] show that a supercritical contact 
process conditioned to survive, when seen on a large scale, stochastically 
dominates a two-dimensional supercritical oriented percolation. They also 
indicate how their construction could be used to obtain a shape theorem. 
This last step essentially consists in proving that the estimates needed 
in |16| hold for the whole supercritical regime, and is done by Durrett [15] 
in 1988. 

Similarly, in the case of a random environment, proving a shape theorem can also 
fall into two different parts. The first one, and undoubtedly the most hard one, 
would be to prove that the growth is of linear order, as soon as survival is possi- 
ble: this corresponds to the Bezuidenhout-Grimmett result in random environment. 
The second one, which we tackle here, is to prove a shape theorem under condi- 
tions assuring that the growth is of linear order: this is the random environment 
analogous of the Durrett-Griffeath work. We thus chose to put conditions on the 
random environment that allow to obtain, with classical techniques, estimates sim- 
ilar to the ones needed in [16] and to focus on the proof of the shape result, which 
already presents serious additional difficulties when compared to the proof in the 
classical case. 

Asymptotic shape results for random growth models are usually proved with the 
theory of subadditive processes, initiated by Hammersley and Welsh [20j, and more 
precisely with Kingman's subadditive ergodic theorem [25] and its extensions. The 
most famous example is the shape result for first passage-percolation on Z d (see also 
different variations of this model: Boivin [5], Garet and Marchand [18], Vahidi-Asl 
and Wierman [35], Howard and Newmann [23], Howard [22], Deijfen |10j). 

The random growth models can be classified in two families. The first and 
most studied one is composed of the permanent models, in which the occupied 
set at time t is non-decreasing and extinction is impossible (such as the Richardson 
models [32], frogs model by Alves and al. [H [2], branching random walks by Comets 
and Popov [9]). In these models, the main part of the work is to prove that the 
growth is of linear order, and the whole convergence result is then obtained by 
subadditivity. 

The second family contains non-permanent models, in which extinction is pos- 
sible. In this case, we rather look for a shape result under conditioning by the 
survival. Hammersley [19] himself, from the beginning of the subadditive theory, 
underlined the difficulties raised by the possibility of extinction. Indeed, if we 
want to prove that the hitting times (t(x)) xeZ d are such that t(nx)/n converges, 
Kingman's theory requires and stationarity and integrability properties for the col- 
lection t(x). Of course, as soon as extinction is possible, the hitting times can be 
infinite. On the other hand, conditioning on the survival can break independence, 
stationarity and even subadditivity properties. A first almost subadditive lemma 
is proposed by Kesten in the discussion of Kingman's paper [25], next improved by 
Hammersley [19] (page 674). Other sets of assumptions are later proposed (see for 
instance Derriennic [TT], Derriennic and Hachem [12], and Schiirger [33l 134] ). 
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The contact process clearly belongs to the second family. Following Bramson 
and Griffeath 0[7], it is on Kesten-Hammersley lemma that Durrett and Griffeath 
rely to prove their shape result. However, their proof contains a certain number 
of inaccuracies that have partially been corrected in [15] . Our strategy, necessarily 
different because of the randomness of the environment, thus offers an alternative 
for the proof of the asymptotic shape theorem for the classical d-dimensional contact 
process. 

Of course, considering a random environment brings extra difficulties. On one 
hand, working in a given - quenched - environment, we loose all the spatial sta- 
tionarity properties. On the other hand, working under the annealed probability, 
we loose the Markovian properties of the contact process. Thus we can not ap- 
ply the Kesten-Hammersley lemma, which requires both stationarity and a kind of 
independence. We introduce here a new quantity <r(x), that can be seen as a regen- 
eration time, and that represents a time when site x is occupied and has infinitely 
many descendants. This a has stationarity and almost subadditive properties that t 
lacks, and that allows to state the problem in the setting of (almost) subadditive 
ergodic theory. We then establish, with techniques inspired from Liggett, a general 
almost subadditive ergodic theorem that gives us an shape result for a. Finally, by 
showing that the gap between t and a is not too large, we transpose to t the shape 
result obtained for a. 



2. Model and results 

2.1. Environment. In the following, we denote by ||.||i and I -lice the norms on M' 
respectively defined by ||x||i = J2t=i \ X A an d IMloo = max \xi\. The notation 

l<i<d 

\\.\\ will be used for an unspecified norm. 

We fix X c (Z d ) < A m in < A max < +oo, where A c (Z d ) stands for the critical param- 
eter for the classical contact process in Z d . In the following, we restrict our study 
to random environments A = (A e ) e6E d taking their value in A = [A m i n , A max ] E . An 
environment is thus a collection A = (A e ) eeE d € A. 

Let A e A be fixed. The contact process (£t)t>o in environment A is a homo- 
geneous Markov process taking its values in the set V(Z d ) of subsets of Z d . For 
x e Z d we also use the random variable £t(z) = l{ 2£ £ t }. If £,t(z) = 1, we say that z 
is occupied or infected, while if £t(z) = 0, we say that z is empty or healthy. The 
evolution of the process is as follows: 

• an occupied site becomes empty at rate 1, 

• an empty site z becomes occupied at rate £t(V)A{ Zj2 '}, 

||z — z'||i=l 

each of these evolutions being independent from the others. In the following, we 
denote by V the set of cadlag functions from R + to "P(Z d ): it is the set of trajectories 
for Markov processes with state space V(Z d ). 

To define the contact process in environment A € A, we use the Harris con- 
struction [21]. It allows to couple contact processes starting from distinct initial 
configurations by building them from a single collection of Poisson measures on R+ . 

2.2. Construction of the Poisson measures. We endow R+ with the Borel a- 
algebra £>(R + ), and we denote by M the set of locally finite counting measures m = 
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S;=o ■ We en dow this set with the er-algebra M generated by the applications 
m i > m(B), where B describes the set of Borel sets in R+. 
We then define the measurable space (f2, T) by setting 

n = M Kd x M zd and T = M m " ® M m " . 

On this space, we consider the family (Pa)asA of probability measures defined as 
follows: for every A = {\ e ) e£ ^d e A, 

Pa= I 0Pa s ) ®V? zd , 

\eGE d ) 

where, for every A £ R+, V\ is the law of a punctual Poisson process on R+ with 
intensity A. If A e R+, we write Pa (rather than P(a) ^ d ) for the law in deterministic 
environment with constant infection rate A. 

For every t > 0, we denote by Tt the a- algebra generated by the applications 
uj h- ► uj e (B) and w i— > uj z (B), where e ranges over all edges in E d , z ranges over all 
points in Z d , and B ranges over the set of Borel set in [0,i]. 

2.3. Graphical construction of the contact process. This construction is ex- 
posed in all details in Harris [H]; we just give here an informal description. Let 

cj = ((w e ) eeE d, (uJ z ) z <£Z, d ) € fi. Above each site z e Z d , we draw a time line R + , 
and we put a cross at the times given by uj z . Above each edge e £ E d , we draw at 
the times given by ui e an horizontal segment between the extremities of the edge. 

An open path follows the time lines above sites - but crossing crosses is forbidden 
- and uses horizontal segments to jump from a time line to a neighboring time line: 
in this description, the evolution of the contact process looks like a percolation 
process, oriented in time but not in space. For x,y e "L d and t > 0, we say that 
£t(v) = 1 if and only if there exists an open path from (x, 0) to (y,t), then we 
define: 

Q - {yeZ d : #(y) = l}, 
(1) VAeV(Z d ) £ t A = |Jtf- 

For instance, we obtain (A C B) (Vi > £f C )• 

When A € R?j_, Harris shows that under Pa, the process (£f)t>o ls the contact 
process with infection rate A, starting from initial configuration A. The proof can 
readily be extended to a non constant A G A, which allows to define the contact 
process in environment A starting from initial configuration A. This is a Feller 
process, and thus benefits from the strong Markov property. 

2.4. Time translations. For t > 0, we define the translation operator 9 t on a 
locally finite counting measure m = J2t=L on R+ by setting 

8 t m = ^ l{u>t}titi-t- 

i=l 

The translation 9 t induces an operator on 0, still denoted by 6 t : for every u e CI, 
we set 

6 t u = ((^tw e ) eeE d, (O t u! z ) ze z d )- 
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The Poisson point process being translation invariant, every probability measure 
Pa is stationary under t . The semi-group property of the contact process has here 
a stronger trajectorial version: for every A c Z d , for every s, t > 0, for every w£S], 
we have 

(2) = d^CM = £(M ° £ t A M, 

that can also be written in the classical markovian way: 

VB e B{V) F((tf +S ) s > G B|^ t ) = P((£).>„ e B) o 

We can write in the same way the strong Markov property: if T is an {Ft)t>a 
stopping time, then, on the event {T < +00}, 

= d #(w) (M, 

VBeBp) P((e# +s ),>o G B\F T ) = P(((;) s >o6B)o? T A 
We recall that .Ft is defined by 

Ft = {B€ J 7 : Vi > Bn{T<t}ef t }. 

2.5. Spatial translations. The group l d can act on the process and on the en- 
vironment. The action on the process changes the observer's point of view of the 
process: for x G Z d , we define the translation operator T x by 

Vo> G ^ T x LO = {(C0x+e) eeK d j K+zlzeZ^i 

where x + e the edge e translated by vector £. 

Besides, we can consider the translated environment x.X defined by (x.A) e = 
\ x +e- These actions are dual in the sense that for every A G A, for every x G Z d , 

(3) Vief V x {T x ueA) = V x . x (LueA). 
Consequently, the law of £ x under Pa coincides with the law of £° under V x .x- 

2.6. Essential hitting times and associated translations. For a set A C Z d , 

we define the life time t a of the process starting from A by 

t a = hd{t > : = 0}- 

For A C Z d and x G Z d , we also define the first infection time t A (x) of site a; from 
set A by 

t A (x) = inf{f >0: 16 

If y G Z d , we write t y (x) instead of t^(x). Similarly, we simply write t(x) for 
t°{x). 

We now introduce the essential hitting time a{x): it is a time where the site x is 
infected from the origin and also has an infinite life time. This essential hitting 
time is defined through a family of stopping times as follows: we set uq(x) = vq(x) = 
and we define recursively two increasing sequences of stopping times (u n (x)) n >o 
and (v n (x)) n >o with uo{x) = vq(x) < Ui(x) < vi(x) < u^{x) . .. as follows: 

• Assume that Vk{x) is defined. We set uu+i{x) = vai{t > Vk(x) : x G ^ }. 
If Vk(x) < +00, then Uk+i{x) is the first time after Vk(x) where site x is 
once again infected; otherwise, Uk+i(x) = +00. 

• Assume that Uk(x) is defined, with k > 1. We set Vk{x) = Uk(x)+r x o6 Uk t x y 
If Uk{x) < +00, The time t x o 9 Uk i x ) is the life time of the contact process 
starting from x at time Uk{x); otherwise, Vk(x) = +00. 
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We then set 

(4) K(x) = min{n > : v n (x) = +oo or u n +x(x) = +00} . 

This quantity represents the number of steps before the success of this process: 
either we stop because we have just found an infinite v n (x), which corresponds to 
a time u n (x) where x is and infected and has infinitely many descendants, or we 
stop because we have just found an infinite u n+ i(x), which says that after v n (x), 
site x is never infected anymore. 

We then set <r{x) = u K r x \, and call it the essential hitting time of x. It is of 
course larger than the hitting time t(x) and can been seen as a regeneration time. 
We will see that K(x) is almost surely finite, so a(x) is well-defined. At the same 
time, we define the operator 9 X on f2 by: 

if <t(x) < +00, 
otherwise, 

or, more explicitly, 

(9 ){uj) = /^(^MM^) if o"(a;)(w) < +00, 
lT x (w) otherwise. 

We will mainly deal with the essential hitting time a(x) that enjoys, unlike t(x), 
some good invariance properties in the survival-conditioned environment. We will 
also control the difference between a{x) et t(x), which will allow us to transpose to 
t(x) the results obtained for a(x). 

2.7. Contact process in the survival-conditioned environment. We now 

have to introduce the random environment. In the following, we fix a probability 
measure v on the sets of environments A = [A m j n , A max ] E . We assume that v is 
stationary and ergodic under the action of X d . This setting naturally contains the 
case of a deterministic environment A > A c (Z d ): we simply take for v the Dirac 
measure (6\)® Ed . 

For A £ A, we define the probability measure Pa on (O, T) by 

MEeT ~P X {E) =P x (E\t° = +00). 

It is thus the law of the family of Poisson point processes, conditioned to the survival 
of the contact process starting from 0. On the same space (fl, T), we defined the 
corresponding annealed probability P by setting 



ME E T P{E) = [ P X (E) dv{\) 

J A 



A 

In other words, the environment A = (A e ) ee i<i where the contact process lives is a 
random variable with law v, and it is under the probability measure P that we seek 
the asymptotic shape theorem. 

It could seem more natural to work with the following probability measure: 

- , N n x [¥x(E)Fx{t° = +oo)dis(\) 

MEET V(E) = F(Er° = +^)= J Al f - >—y . 

It appears that our proofs do not work with this probability measure. However, 
our restrictions on the set A of possible environments ensure that P and P are 
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equivalent: the P— a.s asymptotic shape theorem is thus also aP-a.s asymptotic 
shape theorem. 

2.8. Organization of the paper and results. In Section [3j we establish the 
invariance and ergodicity properties. In particular, we prove the following: 

Theorem 1. For every x £ Z d \{0}, the measure-preserving dynamical system 
(fl,J-,F,9 x ) is ergodic. 

In Section HI we study the integrability properties of the family (a(x)) xe zd ; we 
also control the discrepancy between a(x) and t(x) and the lack of subadditivity 
of a: 

Theorem 2. There exist > such that for any A e A, for any x, y £ Z d , 

(5) Vt>0 W x (<r(x + y)-((r(x)+cr(y)o6 x )>t)<A 5 exp{-B 5 Vi). 

This says that the lack of subadditivity of a is really small: in particular, it does 
not depend on the considered points. Then, in the same spirit as Kingman [26] and 
Liggett [28], we prove in Section [5] that for every x £ Z d , the ratio converges 
P — a.s. to a real number fx(x). The functional x i— > n(x) can be extended into a 
norm on M. d , which will characterize the asymptotic shape In the following, will 
denote the unit ball for We define the sets 

H t = {x£Z d : t(x) < t}, 
G t = {x £ Z d : a{x) < t}, 
K' t = {x£Z d :V S >t e s (x)=tf(x)}, 

and we denote by H t ,G t ,K' t their "fattened" versions: 

H t = H t + [0, G t = G t + [0, l] d and K' t = K' t + [0, l} d . 

We can now state the asymptotic shape result: 

Theorem 3 (Asymptotic shape theorem). For every s > 0, P — a.s., for every t 
large enough, 

(6) (i_ e )^ c _L_* c _* c _* c (l+e)A M . 

The set K' t n Gt is the coupled zone of the process. Usually, the asymptotic 
shape result for the coupled zone is rather expressed in terms of K t C] H t , where 

K t = {x£Z d : £>(*)=£f (*)}. 

Our result also gives the shape theorem for K t C\H t , because K' t C\Gt C K t C\Ht C H t . 

Let us note that the shape result can also be formulated in the following "quenched" 
terms: for v — a.e. environment, we know that on the event "the contact process 
survives", its growth is governed by ([6]) for t large enough. We can also give a 
complete convergence result: 

Theorem 4 (Complete convergence theorem) . For every A £ A, the contact process 
in environment A admits an upper invariant measure m\ defined by 

VA C Z d , \A\ < +oo m x (uj D A) = lim P A (ff D A). 
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Then, for every finite set Actf and for v — a.e. environment X, one has 
Kt Pa(t a < +(x)S +P x (t a = <x)m x , 

where P<J t is the law of £f under P\ and stands for the convergence in law. 

The proof of this result does not require any new idea, and we just give a hint 
at the end of Section [6l 

As explained in the introduction, in order to prove the asymptotic shape theorem, 
we need some estimates analogous to the ones needed in the proof by Durrett and 
Griffeath in the classical case. We set 

B* = {yeZ d : Hi,- *||oo <r} } 

and we write B r instead of . 

Proposition 5. There exist A, B, M, c, p > such that for every A 6 A, for every 
y £ 1 d , for every t > 



(7) P A (r° = +c») > p, 

(8) Px(H?£B M t) < Acxp{-Bt), 

(9) P A (i<T°<+oo) < Aexp(-Bt), 

(10) P\(t°(y)>^- + t, r° = +ooj < Aexpi-Bt), 

(11) P A (0 ^ K' t , t° = +oo) < Aexp(-Bt). 



All these estimates are already available for the classical contact process in the 
supercritical regime. For large A, they are established by Durrett and Griffeath |16j . 
and the extension to the entire supercritical regime is made possible thanks to 
Bezuidenhout and Grimmett's work [4]. For the crucial estimate (9|), one can find 
the detailed proof in Durrett [17] or in Liggett [30]. The need for these estimates 
explains our restrictions on the possible range of the random environment. 

We chose to focus on the stationarity and subadditivity properties of the essential 
hitting time a and on the proof of the shape result. We thus admit in Sections [51 [4] 
and [U the uniform controls given by Proposition [5j whose proof (via restart argu- 
ments) is postponed to Section [6l Section [6] is totally independent of the rest of the 
paper. Finally, in an appendix, we prove a general (almost) subadditive theorem. 
As we think it could also be useful in other situations, we present it in a more 
general form than what is needed for our aim. 

3. Properties of 9 x 

3.1. First properties. We first check that K(x) is almost surely finite and even 
has a sub-geometrical tail: 

Lemma 6. VA c 7L d Vx e Z d VA e A Vn e N F\(K(x) > n) < (1 - p) n . 
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Proof. Remember that p is given in ((7J). Let A G A and n EN. The strong Markov 
property applied at time u n+ i(x) ensures that 

Px(K(x) > n + 1) = F x {u n + 2 (x) < +oo) 

< Px(u n+ i(x) < +oo,v n+ i(x) < +00) 

< Fx(u n+ i(x) < +00, t x o 6 Un+1 ( x ) < +00) 

< P x (u„+i(x) < +cx>)Px(t x < +™) 

< P x (u n+ i(x) < +oo)(l - p) = Px{K{x) > n)(l - p), 
which prove the lemma. □ 
Lemma 7. Let A E A. Pa — a.s., for every x E 1 , 

(12) (K(x) = k) and (r = +00) (uk(x) < +00 and Vk(x) = +00), 

Proof. Let A 6 A. By Lemma [H the number K(x) is Pa — a.s. finite. Let k E N: 
the strong Markov property applied at time Vk(x) ensures that 

Pa(t° = +oo,v k (x) < +oo,u fe+ i(x) = +00 1 F Vk ( x )) 

= ^{v k (x)<+ao}P\{r m =+oo,f(x) = +oo)o£ k{x) . 

Consider now a finite non-empty set B C Z d : with (fT0|l . we get 

P A (r s = +oo,t B (x) = +00) < ^PA(r i/ = +00, ^(x) = +00) 

y£B 

< J2 F y- x ( T ° = +°°,t°(x-y) = = 0. 

yes 

This gives the direct implication. The reverse one comes from ([2]). □ 

Our construction of <j{x) is very similar to the restart process exposed in Durrett- 
Griffeath [16]. The essential difference is that in that paper, the aim is to find, close 
to x, a point that survives while we require here the point to be exactly at x. Thus, 
we will be able to describe precisely the law of the contact process starting from x 
at time a(x), and construct transformations under which P is invariant. 

Lemma 8. Let x E Z d \{0}, A in the a-algebra generated by <j{x), and B E T . 
Then 

VA E A Px{A n (e x Y x {B)) = Px(A)P x .x(B). 

Proof. We just have to check that for any k E N*, one has 

Pa {A n {e x )-\B) n {K{x) = k}) = P X (A n {K(x) = k})P x .x(B). 

Consider a Borel set i'cl such that A = {cr(x) E A'}. The essential hitting time 
a(x) is not a stopping time, but we can use the stopping times of the construction: 

Pa({t° = +00} nAn (e^iB) n {K(x) = k}) 

(13) = P A (r° = +00, a{x) E A', T x o 6 a{x) E B, u k (x) < +00, v k = +00) 

(14) = Px(u k (x) < +00, u k {x) E A', t x o e Uk{x) = +00, T x o 6 Uk{x) e B) 

(15) = Px(u k (x) E A', u k (x) < +oo)Pa(t x = +cx), T x E B) 

(16) = PaK(20 G A', u k (x) < +c^)P x .a({t° = +c^} n B), 

For ifl3|) . we use Equivalence (fl2| . For lfl4|) . we notice that for any stopping time T, 

(17) {T < +00, x E t° o T x o6 T = +00} C {r° = +00}. 



in 
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Equality lfl5|) follows from the strong Markov property applied at time Uk (x) , while 
(HU comes from the spatial translation property (3]). Dividing the identity by 



Pa(t° = +00), we obtain an identity of the following form: 

Pa (A n {ky^B) n {K(x) = k}) = ^{x, A, k, A)F X . X (B), 
and the number tp(x, A, k, A) is identified by taking B = Q. □ 

Corollary 9. Let x, y G 1 d and A € A. Assume that x 7^ 0. 

• The probability measure P is invariant under the translation 9 X . 

• UnderP\, a{y)o9 x anda{x) are independent. Moreover, the law ofa(y)o9 x 
under P\ is the same as the law of o~{y) under F x .\. 

• The variables {cr(x) o (9 x y)j>o are independent under ¥\. 

Proof. For the first point, we just apply the previous lemma with A = £1, then we 
integrate with respect to A and use the stationarity of v. 

For the second point, let A' , B' be two Borel sets in K and apply Lemma [8] with 
A = {a(x) G A'} and B = {a(y) o 9 X G B'}. 

Let n > 1 and An, A±, . . . , A„ be some Borel sets in R. We have: 

P A (a(.T) e A , o6 x &A 1 ,..., a(x) o e A„) 
= P A (a(.T) e A , (a(x), a{x) o o fit, G Ai X • • • X A n ) 

= Pa(ct(.t) G A))Px.a(o-(x) G At,a(x) o 9 X G A 2 , . . . , <r(x) o G A„), 

where the last equality comes from Lemma [U We recursively obtain 

p| Mx) o G A,} = [] P,,.a (a(x) G Aj) , 

,0<j'<« / 0<j<n 

which ends the proof of the lemma. □ 

3.2. Ergodicity. To prove Theorem [IJ it seems natural to estimate the evolution 
with m of the dependence between A and 9~ m (B) for some events A and B. If 
m > 1, the operator 0™ corresponds to a spatial translation by vector ma; and to a 
time translation by vector S m (x): 

9 r x = T mi o 6s m ( x ), 

Tn— 1 

with S m (x) = fi(a;) o 

We begin with a lemma in the same spirit as Lemma [8l 

Lemma 10. Let t > 0, A G .Ft and B E T. 

Then, for any x G Z d , any A G A, any m > 1, 

Pa (A n {t < S m (x)} n {9l n )~\B)) = P A (A n {t < S m (x)})P mx . x (B). 

Proof. Set K m (x) — (K(x),K(x) o 9 X , . . . , K(x) o 9 r x ~ 1 ). Is is sufficient to prove 
that for any k = (k , . . . , fc m _i) G (N*)™ 1 , one has 

Pa (A, t < S m (x), 9 x m (B), K m (x) = k) 
= P A (A, t < S m (x), K m (x) = k)¥ mx . x {B). 
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u k 3 {%) ®l = U k 3 ( x ) ° T 2x ° 0r 2 ( x ) 



U k2 (x) O X = U k2 (x) oT x O Rl ( x ) 



u kl (x) 



x 2x 3x 

Figure 1. An example with ki = 3, k 2 = 2 and k 3 = 4. 



Let k G (N*) m . We set R (x) = and, for I < to - 2, R l+1 (x) = Ri + u kl (x) o 
Or^x). Thanks to Remark (fl7|) . the following sets coincide: 



£° = +00, 

K m (x) = k 



u kl (x) < +oo, Ufc 2 (x) o T x o ^(x) < +oo, 

M fc,„ (x) O T( wi _ 1 ) a! O 0Jj m _ 1 ( a! ) < +oo, 



t u o T_,_ o 6< 



-00 



Moreover, on this event, S m (x) = R m (x) holds. Thus 



r° = +oo, A, 
t < S m (x), 
K m (x) = k, e~ m (B) 



= r\ 



( A, u kl (x) < +oo, 

u k2 (x) oT x o Rl(x) < +oo, 
u k m [x) o T {m _ 1)x o 9 Rm l{x) < +oo, 
i < R m (x), r° o T ma . o 9 Rm(x) = +oo, 



V 



T, 



mi ° "R m (x) 



G B 



J 



By construction, R m (x) is a stopping time and the event 

An {u kl (x) < +00} n • •• n {u km {x) o T( m _ 1)a . o o^^^ < +00} n {i < i2m(x)} 

is measurable with respect to T Rm ^ x y Using the strong Markov property and the 
spatial translation property (J3j) , we get: 



t° = +00, A, 
t < S m (x), 
K m {x)=k, e~ m (B) 



/ A, u kl (x) < +00, \ 

u k2 (x) oT x o 9 Uki (x) < +00, . . . 
u k m {x) oT( ffl -i), o 6' Rm _ 1 ( x ) < +00, 
V t < R m (x) 

xP m , A ({r = +M}nB). 



Dividing the identity by Pa(t = +00), we obtain an identity of the form 

Pa (A t < S m (x), e x m (B), K m (x) = k) = xP(x,\,k,m,A)¥ mx . x (B), 
and we identify the value of ip(x, A, k, to, A) by taking B = Cl. 

We can now state a mixing property: 



□ 
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Lemma 11. Let t > and q > 1 be fixed. There exists a constant A(t, q) such that 
for any x G Z d \{0}, for any A E Tt, for any B G T , A G A and every t>l, 

|P A (An 0i)-\B)) -W x (A)Ft x . x (B)\ < A(t,q)q- e . 

Proof. Let £ > 1. With Lemma flOl we get 

\'P x {Ari0- i (B))-¥ x (A)'Px(^ t (B))\ 
< |F A (t < S e (x), A n - P A (t < ^(z), A)P A (0"*(5))| 

+2P A (i > 
= 2P A (* > Si(x)). 

Let us now fix a > 0. 

With the Markov inequality, P\(Si(x) < t) < exp(ai)E A (exp(— aSe (x))). Using the 
two last points of Corollary [9l one has 

t-i 



E x (exp(-aS e (x))) < E A cxp -a^a{x) o fig 



< JjE A (exp(-a<7(o;) o §1)"} = Y[E jx . x {exp{-aa(x)). 

3=0 3=0 

Now we just have to prove the existence of some a > such that for every A G A, 

E A (exp(— aa(x))) < q^ 1 . 
Let p be the constant given in (7]). 

E A (cxp(— aa(x))) < -E A (exp(— aa[x))) < -E Amsuc (exp(— aa(x))) < 



p p max p a + 2d\ 



max 



because o~{x) > t(x), and t(x) stochastically dominates an exponential random 
variable with parameter 2dA max - This gives the desired inequality if a is large 
enough. □ 

We can now move forward to the proof of the ergodicity properties of the systems 

Proof of Theorem^ We have already seen in Corollary [9] that for any x £ Z d , the 
probability measure P is invariant under the action of 8 X . To prove ergodicity, we use 
an embedding in a larger space to consider simultaneously a random environment 
and a random contact process. 

We thus set = A x fi, equipped with the cr-algebra T = 13(A) <g> J 7 , and we 
define a probability measure Q on T by 

V(A B) G B(A) x T Q(AxB)= [ 1 A (\)P X (B) dv(\). 

J A 

We define the transformation 6a; on Cl by setting 6 X (A, w) = (x.X,8 x (uj)). It is 
easy to see that Q is invariant under Q x . Indeed, for (A, B) G B(A) x T, using 
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Lemma [3 one has 



A 



Q{Q x (\,uj) G A x B) = Q(x.X € A, 9 x (u) € B) 

l A (x.A)P A (flx(w) G B) di/(A) 

1 a (.t.A)P x . a (B) di/(A) 

1 A (A)P A (B) dKA)-Q(AxB) 



Note that if g(X,u) = /(A), then / 3 dQ = / / di/. 
Similarly, if g(A, w) = /(w), then / g dQ = J / dP. 

Note that = U T t is an algebra that generates T . To prove that 9 X is 
ergodic, it is then sufficient to show that for every A e A, 



(18) - £ M^x) converges in L 2 (P) to F(A). 

n k=o 



The quantity above can be seen as a function of the two variables (A, w). Thus, 

n-l 

it is equivalent to prove that the sequence of functions (A, ui) 1— ► i Y] l^f^w) 

fc=0 

converges to P(^4) in L 2 (Q) . Let A G .4 and i > be such that A <E T t . For every 
(u>, A) G f2, we split the sum in two terms: 



1 n— 1 1 n— 1 1 n— 1 

it Z J 77 Z / \ / 77 ^ * 



n ' — ' n * — ' v / n 

k=0 k=0 k=0 



If we set /(A) = Pa (A), the second term can be written 

n— 1 n—1 

-VUA) = T/(M. 

77 Z J nri ' J 



n i — ' n 

k=o k=o 



Since v is ergodic, the Von Neumann ergodic theorem says that this quantity con- 
verges in L 2 {v) to J fdv = ¥(A). Seen as a function of (A, w), it also converges in 
L 2 (Q) to f(A). Set, k > 0, 
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and L n = Y + Y\ + ■ ■ ■ + Y n -\. It only remains to prove that L n /n converges to 
in L 2 (Q). As Ffe = Y o 0^, the field {Y k ) k > a is stationary. We thus have 



L 



< 



< 



E 

0<i,j<n— 1 
n-1 



YiYj 



IL J- n It, ± „ 

]T \ Yi dQ + 2 Y,( n - £ ) / YoYe dQ 

a n *f d i " 

(+oo - N 
J2 / \Kx(Y a Y e )\ du(X) 

1-*(A)) -¥ x (A)¥ x (e-«(A))\ du(X)\ 



(+oo . 
^ / 
1=0 J 

(+oo „ 
^ L 
e=o JA 



Y Y e 



' +oo 



< 2n [Y A(t, 2)2- e = 4A(t, 2)n, 



\i=0 



thanks to Lemma[TTJ This ends the proof of (fl8)) , hence the proof of Theorem[TJ □ 



4. Control of the subadditivity flaw 

In this section, we are going to bound quantities such as a(x+y) — [a(x)+a(y)o9 x ] 
and <j{x) — t(x). 

We will use these results in the application of a (almost) subadditive ergodic 
theorem in Section[5l In both case, we use a kind of restart argument. Considering 
the definition of the essential hitting time a, we will have to deal with two types 
of sums of random variables, that are quite different: sums of Uj — Ui on one hand, 
and sums of — Vi on the other hand: 

• The life time Vi{x) — Ui{x) of the contact process starting from x at time 
Ui(x) can be bounded independently of the precise configuration of the 
process at time Ui(x), So the control is quite simple. 

• On the contrary, — Vi(x), which represents the amount of time 
needed to reinfect site x after time Vi(x), clearly depends on the whole 
configuration of the process at time Vi(x), which is not easy to control 
precisely and uniformly in x. This explains why the restart argument we 
use is more complex and the estimates we obtain less accurate than in more 
classical situations (For instance, in Section [6l we obtain the exponential 
estimates of Proposition [5] by standard restart arguments.) 

As an illustration of the first point, we easily obtain: 

Lemma 12. There exist A,B>0 such that for every X G A, 

(19) VxeZ d Vi>0 P\(3i < K(x) : v l (x) - u,(x) > t) < Aexp(-Bt). 

Proof. Let F : — > R + be a measurable function and x E Z d . We set 

£-x(F) = l{ Ui (x)<+oo}^ ° Ouj(x)- 

i=0 
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With the Markov property and the definition of K (x), we have 

+00 / +00 \ 

E X [C X (F)} = ^E A [l K(:c )< +oo} ]E A [F]= ll + J2V\(K(x)>i) \E X [F] 

i=0 \ i=0 / 

= (1 + E x [K(x)])E x [F] < (l + -^j E X [F], 

where the last equality comes from LemmalU We choose F = TL{t<ui(x)-vi(x)<+oo}, 
and with Estimate ([7]), we obtain: 

P x (3i < K(x) : vAx) - uAx) > t) < -F x (3i < K(x) : vAx) - mix) > t) 

P 

< -P X (C X (F) > 1) < -E X [C X (F)) 
P P 

< ~ ( 1 + M P A (t < T X < +00). 

We can then conclude with Inequality ([9]). □ 

To deal with the reinfection times Ui + \(x) — Vi(x), the idea is to look for a 
point (y, t) - in space-time coordinates - close to (x, Ui(x)), infected from (0, 0) and 
with infinite life time: The at-least-linear-growth estimate iflOf will then ensure it 
does not take too long to reinfect x after time t, just by looking infection starting 
from the new source point {y,t). The difficulty lies in the control of the distance 
between (x,Ui(x)) and a source point (y,i): if the configuration around (x,Ui(x)) is 
"reasonable", this point will not be too not far from (x,Ui(x)), and we will obtain 
a good control of Ui + i(x) and Ui(x). 

We recall that for every x G Z , u x is the Poisson point process giving the 
possible death times at site x, and that M and c are respectively given in ([8]) 
et (TlOl). Note that we can assume that M > 1. We note 



(20) 7 = 3M(1 + 1/c) > 3. 

For x, y € Z d and t > 0, we say that the growth from (y, 0) is bad at scale t with 
respect to x if the following event occur: 

E»(x,t) = K[0,i/2] =0}U{ff t y £y + B Mt } 

U {t/2 <t v < +00} U {t v = +00, inf{s > 2t : i€ ^} > 7*}, 

We want to check that with a high probability, there is no such bad growth point 
in a box around x: for every x £ Z d , every L > and every t > 0. So we define 

In other words, we count the number of points (y, s) in the space-time box {x + 
Bah+i) x [0, L] such that there happens something for site y at time t, either a 
possible death, or a possible infection, and at this time the bad event E y {x,t) o 6 S 
occurs. We first check that if the space-time box has no bad points and if Ui(x) is 
in the time window, then we can control the delay before the next infection : 

Lemma 13. // Nl(x, t) o 9 S = and s + t < Ui(x) < s + L, then Vi(x) = +00 or 
u i+1 {x) - Ui(x) < -ft. 
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Proof. By definition of Ui(x), site x is infected from (0,0) at time Ui(x). Since 
s + t < Ui(x) < s + L and Ui(x) is a possible infection time for x, the non-occurrence 
of E x {x,t) o 9 Ui (~\ ensures that r x o 9 Ui r x \ = +00 or that r x o 9 Ui r x \ < t/2. If 
t x o 6 Ui ( x ) = +00, we are done because then Vi(x) = +00. Otherwise, note that 
Vi(x) — Ui(x) < t/2. 

By definition, there exists an infection path 7, : [0, Ui(x)] — > Z d from (0,0) to 
{x,Ui(x)), i.e. such that 7i(0) = and Ji(ui(x)) = x. Consider the portion of 7, 
between time Ui(x) — t and time Ui(x). Denote by x = ^n(ui(x) — t) and let us see 
that xo £ x + Bnjt+2- Indeed, if xo </ x + Bmt+2, we seek the first time ti after 
time Ui(x) — t when 7* enters in x + B^t+2 at a site we call x\ (note that since x\ 
is in the inside boundary of x + Bj^t+2, we have ||x — XiHqo > Mt + 1): time t\ is 
a possible infection time for Xi, and the non-occurrence of E Xl {x,t) o 9 tl ensures 
that the infection of x from (xi, ti) will at least require a delay t, which contradicts 
Ui(x) -t>0. 

So xo G x + BMt+2- since Nz(x,t) o 9 S = 0, the first possible death at site xo 
after time itj(x) — t can not occur after a delay of t/2; thus the first time £2 when 
the path 7, jumps to a different point X2 satisfies ta < u-s(x) — t + t/2 = Ui(x) —t/2. 
Consequently, when {x2,t 2 ) infects (x,iti(x)), it is at least t/2 aged, and the non- 
occurrence of E X2 (x, t) o 9 t2 ensure it lives forever and 

inf{w>2t: xeC 2 }o9 t2 < 7 t. 

So there exists t 3 € [t 2 + 2t,t 2 + 7*], with x G £° 3 . Since V,(x) — u,(x) < t/2, one 
has 

t3 > t 2 + 2t > (u;(x) - t) + 2t = m(x) + t > Vi{x). 
Finally, Uj_|_i(x) — Uj(x) < t% — u%(x) < t2 — Ui{x) + jt < 7t. □ 

Now we estimate the probability that a space-time box contains no bad points: 

Lemma 14. There exist Aref^ ^2l\ ■ > such that for every A € A, 

(21) VL>0 VxeZ d Vt>0 P A (JVi(x,t) > 1) < A 21 (l + L)exp(-B 21 t). 
Proof. Let us first prove there exist A, B > such that for every A G A, 

(22) Vx e Z d Vt > Vy G x + S A/t+2 P A (i7 y (x,t)) < Aexp(-St). 

Let x G Z d , t > and y G x + BMt+2- If r 2 ' = +00, there exists z G £ft with 
r z o 6>2* = +00. Thus the definition (|20| of 7 implies that 

{r y = +00, inf{s > 2t : x G 0} > 7^} 

C {&2y + fl 2 Jift}U |J {t z (x)o0 2t > ( 7 -2A/)*} 

z£y+B2Mt 

c UI^ y + ^MjU (J lt z (x)o9 2t > h^A\ +M t-- c 

z£y+B 2 Mt 

Hence, with © and ((10)1 . 

Pa(t !/ = +00, inf{s > 2t : x G £*} > 7*) 

< Acxp(-2BMt) + (1 + 4Aft) d Acxp(-.B(Mt - 3/c)). 

The distribution of the number w y ([0,t/2]) of possible deaths on site y between 
time and time t/2 is a Poisson law with parameter t/2, so 

P A H([0,t/2]) = 0) = exp(-t/2). 
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The two remaining terms are controlled with ([8]) and ([9]); this gives ([22]) . 

Fix now y G x + and note f3 y = tu y + lu £ . Under Pa, f3 y is a 

Poisson point process with intensity 2d\ e . Let So = and (S n ) n >i be the increasing 
sequence of the times given by this process. 

f L +°° 

/ ^Ey(x.t) ° &s d(P y + Sq)(s) = ^ l{S„<L}l£»(x,t) ° &S„ ; 

so, with the Markov property, 

-j-oo -\-oo 

= ]T E A (l {Sn <L}lEv( Xl t) O 9 Sn ]) = J2 E A (1{5„<L } ) Pa(^(^, 0) 
n=0 ra=0 

= (l + E A [^([0,L])])PA(S' y (.T,<)) = (l + L(2d\ e + l))P x (Ey(x,t)). 

So <(2TJ follows from ([22]), from the remark that Vx(N L (x,t) > 1) < E A [iV L (a;, t)], 
and from on obvious bound on the cardinality of But+2- □ 

Once the process is initiated, Lemma [13] can be used recursively to control 
Ui + \(x) — Ui(x). To initiate the process, we assume that there exists a point (u, s), 
reached from (0,0), living infinitely and close to x in space: 

Lemma 15. For any t, s > 0, for every x G Z d , i/ze following inclusion holds: 

{t° = +00} n{3uex + B Mt+2 , t u o6 s = +00, u g £° } 

(23) n{iV ir(;E)74 (a ; ,t)o^=0} 

(24) n pi {vi{x) - Ui (x) < t} 

l<i<K(x) 

(25) C {t° = +00} n {a(x) <s + K(x)yt} . 

Proof. If every finite Ui{x) is smaller than s + t, we are done because a(x) < s + t < 
s + K(x)yt. So set 

io = max{i : Ui(x) < s + t}. 

Since Vi (x) < +00, the event ([24]) ensures that «i (a;) — Ui (x) < i, and so Vi (x) < 
s + 2t. Now, since r" = +00, the non-occurrence of E u {x,t) o 8 S implied by ([23]) 
says that 

inf{s > 2i : x G £"} ofl s < -yt, 

which leads to itj +i(x) < s + 'yt. Noting that for any j > 1, Ui 0+ j(a- + y) > s + t, 
we prove by a recursive use of Lemma fT3l with the event {N K ^ lt (x,t) o 9 S = 0} 
that 

Vj G {1, . . . , X(x) - i } u io+ j < s + i7*. 

For j = K(x) — io, we get o~{x) = Ui 0+ j(x) < s + (K(x) — ^0)7^ < s + K(x)jt, which 
proves ([25]). □ 
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4.1. Bound for the lack of subadditivity. To bound a[x+y) — \<y{x)+(x{y)od x \, 
we apply the strategy we have just explained around site x + y. To initiate the 
recursive process, one can benefit here from the existence of an infinite start at the 
precise point (x + y, <j(x) + a(y) o 9 y ). 

Proof of Theorem® Let x, y E Z d , A E A and t > 0. We set s = a(x) + a(y) o 6 X . 
Px(<r(x + y)> a(x) + <r(y) o § x + t) 

Vt\^J t° = +oc, K(x + y)<f \ 



< P A K{x + y)> 



7 / \ <j{x + y) > s + K(x + y)7v / i 



With the sub-geometrical behavior of the tail of K given in Lemma [6] and the 
uniform control ([7]), we can control the first term. Note that if K{x + y) < 
then K(x + y)jVt < t, and so that 

{N K{x+y)lVI (x + y,Vi)>l}c{N t (x + y,Vt) > 1}. 

We apply Lemma [15] around x + y, on a scale \/t, an initial time s = a(x) + a(y)o9 x 
and a source point u = x + y: 

t° = +oo, K(x + y) < — , a{x + y) > s + K(x + y)"/VtJ 

< Px(N t (x + y,Vt)o9 s >l) 

(26) +V x (3i < K(x + y) : v l {x + y) - m{x + y)>\Tt). 

Since N t (x + y, ^/t) = N t (0, Vt) o T x o T y and s = a(x) + <j(y) o 8 X , we have 

N t (x + y,V~t)o9 s = N t (0, Vt)o6 y o 9 X . 

Thus W\(N t (x + y,Vi)o 6 S > 1) = ¥( x+y ). x (Nt(0, Vi) > 1), which is controlled by 
Lemma fl4l and estimate ([7]). Finally, |26|) is bounded with Lemma fl2l □ 

Corollary 16. For x, y E Z d , set r(x, y) = (cr(x + y) — {<j(x) + <j(y) o 9 X )) + . 
For any p > 1, there exists M p > such that 

(27) VA E A Vx, y E 7L d E x [r{x, yf] < M p . 

Proof. We write E\[r(x, y) p ] — J^°° pu p ~~ 1 F\(r(x,y) > u) du and use Theorem^ 

□ 

4.2. Control of the discrepancy between hitting times and essential hit- 
ting times. To bound a(x)—t(x), we would like to apply the same strategy starting 
from (x,t(x)) but we do not have any natural candidate for an infinite start close 
to this point. We are going to look for such a point along the infection path be- 
tween (0,0) and (x,t(x)), which requires controls on a space-time box whose height 
- in time - of order t(x), i.e of order ||a;||. So we will lose in the precision of the 
estimates and in their uniformity. 

Proposition 17. There exist ^gg) , c fJ2l ^ ^ such that for every z > 0, every 
x E 7L d , every A E A; 

(28) P A (<j(x) > t(x) + K(x)(a 2S log(l + ||x||) + z)) < A 28 cxp(- J B 28 z). 
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Proof. For x,y £ 7L d and t, L > 0, we define 

E«(t) = {t v < +00, L) n H x s gL B Mt }, 



s>0 



N L (x,t) = 

yex+B M t+ 



/ *-EV{t)° B » d \ E UJe ) 



With ([7]), © and (({j}, it is easy to get the existence of A, B > such that 

(29) VA £ A V.t e Z d V<> ¥ x (N L {x,t) > 1) < A(l + L) cxp(-Bt). 

Now, we choose the last point (u, s) on the infection path between (0, 0) and {x, t(x)) 
such that t u o 9 S = +oo. Note that on {r° = +oo}, such an s always exists. 

Let us see that if N t ^(x, t) = 0, then u £ x+B Mt+2 . Indeed, if x\\ > AIt+2, 
we consider the first point (u', s') on the infection path after (u, s) to be in x + BMt- 
The definition of s ensures that the contact process starting from (u' , s') does not 
survive, but, since it contains (x,t(x)), its diameter must be larger than Mt, which 
implies that N t ( x ){x,t) > 1, and gives the desired implication. 

On event {N t ( x )(x, t) = 0}, we are going to apply Lemma fl5l around point (x, 0), 
at scale 

= alog(l + ||x||) + z > z 
1 ~ 7' 

with source point («, s) and a time length L — K(x) r yt. Here and in the following, 
a > is a large constant that will be chosen later. Since s < t(x), 

Pa (<j(x) > t{x) + K(x)(a\og(l + \\x\\) + z)) = P A (a(x) > t{x) + K(x)"/t) 

< Pa (a(x) >s + K(x)jt) 

< Pa (a{x) > s + K(x)jt, N t(x) (x,t)=0j + P A (N t{x) (x, t) > l) 

(30£ F x (N K(x)lt (x,t) o 6 S > 1) + P x {3i < K(x) : Vi (x) - Ui (x) > t) 
+¥ x (N t{x) (x,t) > 1). 

The second term in l|30p is bounded with Lemma fT2l For the last term, we write: 

¥x(Nt(x)(x,t)>l) < P x (N^ i+z (x > t)>l)+¥x(t(x)>^ + zj . 
The second term is controlled with ([7]) and (jTUfl . and l|29p ensures that 
Pa (#j[£ji +t (ar, t) > l) < A f 1 + M + ^ cxp(-Bi) 

< A (l + M +Z ] exp f B(alo g (l + M ) +Z ) 



c J \ 7 
< A'exp(-B'z), 

as soon as a is large enough. 

For the first term of |30|) . we note that N K ^ lt (x, t)o6 s < N t ^ +K ^ lt {x, t). Thus 
¥x(N K{xht (x,t)o6 s > 1) < Pa (N^ +z+K(xht (x,t) > l) + P A (t(x) > M + zj . 
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As previously, the second term is bounded with ([7]) and (fTO)) . while using l|2Tj) . we 
get: 

Pa (N^ +Z+K{xht (x,t)>l 



fc=l 
+ OO 



^ \ AlH f 1 + ^r) (! + ^ + ^) ex P f £ V(l + fc)PA(A-(a:) = fc). 



fc=i 



The sub-geometrical behavior of the tail of K(x) given by Lemma ([6]) ensures that 
the sum is finite, and we end the proof by increasing a if necessary. □ 

Lemma 18. For every p > 1, there exists C^j^p) > such that for every x e Z d 

(31) VAeA E x (\a(x)-t(x)n<C 31 (p)(log(l + \\x\\)) p . 

Proof. Set V x = ^fepf^ — of28ll°g(l + IMD- By Proposition [T7l there exists a 



K(x) 

random variable W with exponential moments that stochastically dominates V x 
under Pa for every x and every A. Moreover, Lemma [6] ensures that K(x) is 
stochastically dominated by a geometrical random variable K' . 

Set v(x) = a(x) - t(x) = K(x)(alog(l + ||x||) + V x ) and let p > 1. With the 
Minkowski inequality, we have 

(E x v{x) p ) 1/P < a\og(l + \\x\\)(E x K(x) p ) 1 / p + (E x [K(x) p VV]) 1/p 

< alog(l + \\x\\)(E x K(x) p )^ p + (E x K(x) 2p E x V x 2p )^ 

< alog(l + ||a;||)(Ei<" p ) 1/p + (EK /2p EW 2p ) , 

which ends the proof. □ 

Corollary 19. P - a.s., hm ^ x ]~ ^ = . 

||z||— »+oo \\.r\\ 

Proof. Let p > d: Equation lj3Tj) gives 

w Kx)-twr ^ (iog(i + iixi))" 

So ( ^(^[jjf)^ J d i s almost surely in £ p (Z d ), and thus goes to 0. □ 

Corollary 20. There exist ^gg], -E fggt ^tsE ^ ^ SMC ^ / or every A G A, 
(32) Vx e Z d Vi > P A (cr(x) > C32 1|^|| + t) < ^32 exp(-B 3 2\/t). 



Proof. Let a = a|2g]be given in Proposition[17j and note that if A"(x) < ^ 



and 2 = a;^/||x|| + t/2, then, since log(l + u) < ^Ju, we get: 

/if(x)[olog(H-||s||) + z] <2zK(x)< \\x\\+t/2. 



ASYMPTOTIC SHAPE FOR THE CONTACT PROCESS IN RANDOM ENVIRONMENT 21 



Thus with HH and Q, 

Pa (^0*0 > Q + 1 1 ^1 1 ^ * 
< Pa > M + t/2) + Pa ( (x) > ^VW\+tJ^ 

+P A (cr(a-) > t(x) + K(x){alog(l + \\x\\) + ay/\\x\\ + t/2)} . 

The first term is controlled with lflO|) . the second one with Lemma [6j and the last 
one by Proposition [17l □ 

Corollary 21. For any p > 1, there exists £fgg|(p) > such that 

(33) VAeA V.TGZ d E x [a(x) p }<C 33 {p)(l + \\x\\) p . 

Proof. With the Minkowski inequality, one has 

(IaH*) p ]) 1/p < q^lMI + (Ea[((^(x) - q^MI) + n) 1/p - 

Moreover 

/■+00 

Ea[((<t(x) - C^\x\\)+r] = / ^- 1 P A (a(x) - C^\x\\ >u)du< +^ 
by Corollary EQl □ 



Remarks. In classical restart arguments, the existence of exponential moments 
for a random variable usually comes from the following argument: if (X„)„ 6 pj are 
independent identically distributed random variables with exponential moments, 
if K is independent of the {X n ) n ^s, and also has exponential moments, then 

X n has exponential moments. Here, our difficulties to precisely bound the 

0<n<K 

reinfection times Uj+i — Vi prevents us to use this scheme: we thus have to use ad 
hoc arguments, which lead to weaker estimates. 

5. Asymptotic shape theorems 

We can now move forward to the proof of Theorem [31 The first step consists in 
proving convergence for ratios of the type "'(""v , With Corollary [TBI we know that 
for every n,p> 0: 

E[a((n+p)x)} < E[a(nx)] + E[a(px)} + Mi. 

Thus the Fekete lemma says that iE[cr(nx)] has a finite limit when n goes to +oo 
and the natural candidate for the limit of SlltEl [ s thus 

n 

E(cr(nx)) 
fj,(x) = hm . 

n—*+oo n 

™ ™rf a(nx) , Eer(rax) , . 
Theorem 22. P - a.s. Vx e 7L d lim — — - = lim ^ — '- = fi(x) . 

This convergence also holds in any L P (P), p > 1. 

To prove this result, we need the two following (almost) subadditive ergodic 
theorems, whose proof will be given in the appendix. 
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Theorem 23. Let (f2, .F, P) be a probability space, {9 n ) n >i a collection of trans- 
formations leaving the probability measure P invariant. On this space, we consider 
a collection (f n )n>i of integrable functions, a collection (g n ) n >i of nonnegative 
functions, and a collection (r n . p ) n . p >i of real functions such that 

(34) \/n,p>l f n+P <f n + f p o 

We assume that 

• c = inf„>i ^ > -oo. 

• gi is integrable, g n /n almost surely converges to and converges to 0. 

• There exists a > 1 and a sequence of positive numbers {C p ) p >\ such that 
^[( r «.p) Q ] — f or every n,p and 

+oo n 

V^<+oo. 
Then ±E/„ converges; if p denotes its limit, one has 



E 



1 ■ fn 



n — >+oo 71 

If we set f = lim then f is invariant under the action of each 6 n . 

n — >+oo 

Theorem 24. We keep the setting and assumptions of Theorem \23[ We assume 
moreover that for every k, 

^fnk-Yjk°(O k ) 1 ) ^0 a.s. 
Then f n /n converges a.s. to f. 



i=0 



Proof of Theorem\2^ We apply Theorem[23]with the choices /„ = a(nx), 6 n = 9 nx , 
g p = 0, r n:P = r(nx,px) and the probability measure P = P. We take a > 1. 
Corollary [21] gives the integrability of a(x) under P and Corollary [16] gives the 
necessary controls on its moments. 

We now check the extra assumption of Theorem [24j it is easy to see that 



n-l 



t(nkx) < a(kx) o (6 kx ) 1 



i=0 



n-l 



which implies that I a(nkx) — a(kx) o {6~k x ) 1 ) < a(nkx) — t(nkx); 

V t=o / 

Corollarv[19l ensures that this quantity is o(n). Thus a(nx)/n converges to a random 

variable fj.(x), which is invariant under the action of 6 X . But Theorem [T] says that 

this fx(x) is in fact a constant, which ends the proof of the a.s. convergence. 

To prove that a sequence converges in L p , it suffices to show that it converges 

a.s. and that it is bounded in L q for some q > p. Since Corollary [21] says that f n /n 

is bounded in any LP , the proof is over. □ 

The next step is to prove the asymptotic shape result, namely Theorem [3] We 
start by proving the shape result for the essential hitting time a, by following the 
classical strategy: 
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• We extend /i to a norm on R d in Lemma l25l 

• We prove that the directional convergence given by Theorem [22] is in fact 
uniform in the direction in Lemma [27l 

• We easily deduce the shape result from this lemma in Lemma l28l 

To transpose this shape result for the classical hitting time t (Lemma [29]), we 
just need to control the discrepancy between a and t: this was done in Lemma [191 
Finally, the shape result for the coupled zone is proved in Lemma[30l by introducing 
a coupling time t' and by bounding the difference between this time t' and the 
essential hitting time a. 

Lemma 25. The functional fi can be extended to a norm on M. d . 

Proof. Homogeneity in integers. We know that fi(x) = limE °^ ng ^ , and that a(nx) 
and a(—nx) have the same law under P, so fi(x) = u(—x). By extracting subse- 
quences, we prove the homogeneity in integers: 

Vfc e z Vx e z d n{kx) = \k\fj,(x). 

Subadditivity. One has a(nx + ny) < a(nx) + a(ny) o 9 nx + r(nx, ny). 
Since P is invariant under the action of nx , we get with Corollary [TBI 

Ea(na; + ny) < ¥,o~(nx) + Ea(ny) + Er(nx, ny) < ~Ea(nx) + Ecr(ny) + Mi. 

We deduce that Va; E Z d My e Z d u(x + y) < fi(x) + n(y). 

Extension to R d . The Fekete Lemma ensures that 

, , . , Ea(nx) + Mi 

a(x)+M 1 = mf — -, 

«>i n 

so u(x) < Ea(x). Corollary [21] gives some L > such that Ea(x) < L\\x\\ for any 
x. Finally, n(x) < L\\x\\ for every x £ Z d , which leads to \n(x) — fJ.{y)\ < L\\x — y\\ : 
we can then extend n to Q d par homogeneity, then to R d by uniform continuity. 

Positivity. Let M be given by Proposition [5] With |(8]), we obtain 

~ J P A (r° = +oo) 

A ( n\\x\\ 
< — cxp —B 



p \ 2M 

With the Borel-Cantelli Lemma, we deduce that fi(x) > ^jj\\x\\. This inequality, 
once established for every x £ Z d , can be extended by homogeneity and continuity 
to R d . So n is a norm. □ 

In the following, we set C = 2Cj^2j where Cj^in given in Corollary [20l 
Lemma 26. For every e > 0, P — a.s., there exists R > such that 

Vx,y e 7L d (||a;|| > i? and \\x-y\\ < e\\x\\) (\a(x) - a(y)\ < Ce\\x\\). 
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Proof. For m G N and e > 0, we define the event 

A m (e) = {3x, y £ 7L d : \\x\\ = m, \\x — y\\ < em et \cr(x) — a(y)\ > Cem] 
Noting that 

A m {e) C (J {a(y-x)o6 x +r(x,y-x)>Csmj, 

(l-e)m<||a;||<(l+e)m 
||a; — y||<em 

we see, with Corollaries [20] and [T6] that 



F x (A m (e)) < 



< 



^2 p A (<r(z) °0x + r(x, z) > Cem) 
=IMI< 

Z ||< E 7 

]T P x .a(<t(2) > 2Cem/3) 



(l-e)m<||a:||<(l+e)m 
|k||<em 



(1 — e)m<||a:||<(l+e)m 
|U||<em 



< 



+P x (r(x,y - x) > Cem/3) 
(1 + 2em) d (l + 2(1 + e )m) d ^ exp(-E^/Cem/3) 
+^cxp(-£j27]v / CW3)) 

by Corollary [20] and Theorem [2] Integrating then with respect to A, we conclude 
the proof with the Borel-Cantelli Lemma. □ 

— , \cr(x) — u(x)\ 

Lemma 27. P - a.s. Urn 1 v ; „ „ v n = 0. 

11*11—+°° ||x|| 

Proof. Assume by contradition that there exists e > such that the event "|cr(ar) — 
fi{x)\ > e\\x\\ for infinitely many values of x" has a positive probability. We focus on 
this event. There exists a random sequence (y n )n>o of sites in 1 d such that \\y n \\i — > 
+oo and, for every n, |cr(y n ) — /i(y n )| > e||y rl |ji. By extracting a subsequence, we 
can assume that 

WVnh ' 

Fix £\ > (to be chosen later); we can find z' £ 1 d such that 

z' 

< ei. 



\z'\\i 

For each y n , we can find an integer point on Rz' close to y n . Let h n be the integer 
part of Iffiii ^ ; we have 



\\y n — h„.z'\\i < 

< llVnlll 



||lfa||l ; 



llj/n 



z 1 



y« i ll z Hi 



Take AT > large enough to have (n > iVj =>■ (| 
for z' , one has 



: 'lll 
+ \Wh 

— z \\i — £ i)- By our choice 



(n > N) 



\Vn\\l \\Z'\\l 



< 2ei , 
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and, consequently, \\y n — h n .z'\\i < 2ei\\y n \\i + \\z'\\i. Thus, increasing N if neces- 
sary, one has, for every n> N, \\y n — h n .z'\\ i < 3ei||y„||i. But if AT is large enough, 
Lemma l26l ensures that 

Vn > N \a{y n ) - a{h n .z')\ < ZCe x \\y n \\x. 

Finally, for every large n we have 

WiVn) - ^{Vn)\ < W(Vn) - cr(h n .z')\ + \a(h n .z') - fi(h n .z')\ + \n{hn.z') - n{y n )\ 

a{h n .z') 



< 3Csi\\y n \\i + h Ti 



< 3Cei||y„||i + (l + ei 



h„ 

huh 



- n(z') 



L\\h n .z' - y n \\ x 



a(h n .z') 



H(z'] 



3e-iL||y„ 



But the a.s. convergence in the z 1 direction ensures that for every large n, 

a(h n .z') 



H(z') 



< 



Now if E\ > is small, we obtain, for every large n, \a(y n ) — fi(y n )\ < e||2/n||i; this 
brings contradiction and ends the proof. □ 

We can now prove the shape result for the "fattened" version Gt of Gt — {x G 
Z d : <j{x) < t}; we recall that is the unit ball for fx. 

Lemma 28. For every e > 0, P — a.s., for every large t, 

0—e)A„c ^c(l+e)A^ 

Proof. Let us prove by contradiction that if t is large enough, -j- C (1 + s)A )Ji . 

Thus assume that there exists a increasing sequence(i n )„>i, with t n 

G t 



and 



-p 4 - <£_ (1 + e)A ll : so there exists x„ with a(x n ) < t n and n(x n )/t n > 1 + e. So 
n(x n )/a(x n ) > 1 + £, which contradicts the uniform convergence of Lemma l27l 
since fi(x n ) > t n (l +e), the sequence (||x„||) n >i goes to infinity. 

For the inverse inclusion, we still assume by contradiction that there exists a 

increasing sequence (< n )n>i, with t n — > +oo and (1 — e)A^ <f_ -f^] this means we 
can find x n with fi(x n ) < (1 — s)t n , but a(x n ) > t n . Since t n goes to +00, the 
sequence (x n )n>i is not bounded and satisfies ^ Xn \ < 1 — e; this contradicts once 
again the uniform convergence of Lemma ED and ends the proof. □ 



Then we immediately recover the uniform convergence result for the hitting 
time t via Lemmafl9l and, by an argument similar to the one used in LemmaEHl the 
asymptotic shape result for the "fattened" version Ht of H t = {x E 7L d : t(x) < t}. 



Lemma 29. 



and for every s > 0, 



lim 

II — ~ + °° 



t(x) — fx(x) 



0, 



Hi 



a.s., for every large t, (1 — e)A fl C C (1 + e)A fl . 



It only remains now to prove the shape result for the coupled zone K[, which is 
the "fattened" version of K[ = {x £ Z d : Ws>t = ^(a;)}: 

K r n G 

Lemma 30. For every e > 0, P — a.s., for every large t, (1 — e)A fl C — — . 
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Proof. Since 1 1— > K' t fl G t is non-decreasing, we use the same scheme of proof as for 
Lemma EH We set, for x E Z d , 

t'{x) = m£{t > : x E K' t H G t }. 
It is then sufficient to prove that P — a.s., lim |t (x)-a(x)\ _ q 

\\x\\->+oo 

By definition, t'(x) > a(x); thus it is sufficient to prove the existence of constants 
A', B' > such that 



(35) Vx G 1 d Ms > ¥(t'(x) - a(x) > s) < A'e 



-B's 



• Note first that for every t > 0, K^m^ x + K t o 9 x . 

Indeed, let z E x + K t o 8 X . Consider first the case z £ ^r x )+f Since, by additiv- 

ity ®, C(x)+t c CU+t> we have 2 £ C(x)+f and so that z € ^(x)+f 

Consider now the case z E & x \ +t - Since, by additivity, & x -\ C £o ° we 
have y = z — x E £f o 9 X . But since y € K t o X , the definition of -Kt implies 
that £ t °(y) o X = £? d (y) o ^ = 1. Since x G £° (x) and y G ^ o X , we obtain 
z = x + y E C(x)+v and so z e K <y(x)+t- 

• Fix s > 0. The previous point says that 

f) K a(x)+t J D I x + f\(K t o 61.) , and so D (x + o a )) . 

Since P is invariant under 8 X , we get: 

¥(t'(x)>a(x) + s) = P(x i K' a{x)+S n G a(:c)+S ) 

= P(:r g /C,, )+s ) 

< ¥(x(£(x + K' s o6 x f) 

< ¥{0(tK' s o9 x )=¥(0(£K' s ). 

We conclude with (11]) • □ 



6. Uniform controls of the growth 

The aim of this section is to establish some of the uniform controls announced in 
Proposition [H To control the growth of the contact process, we need some lemmas 
on the Richardson model. 

6.1. Some lemmas on the Richardson model. We call Richardon model with 
parameter A the time-homogeneous, ^(Z^-valued Markov process (r]t)t>o, whose 
evolution is defined as follows: an empty site z becomes infected at rate A "V^ r\t (z'), 

\\z-z'\\ 1 = l 

the different evolutions being independent. Thanks to the graphical construction, 
we can, for each A E A, build a coupling of the contact process in environment A 
with the Richardson model with parameter A ma x, in the following way: at any time 
t, the space occupied by the contact process is contained in the space occupied by 
the Richardson model. 
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The first lemma, whose proof is omitted, easily follows from the representation 
of the Richardson model in terms of first passage percolation, together with a path 
counting argument. 

Lemma 31. For every A > 0, there exist constants A, B > such that 

Vt>0 P(j7i £ B t ) < Aexp(-Bt). 

Lemma 32. For every A > 0, there exist constant A, B, M > such that 

Vs > P(3f > : 7? t (fL B M t+s) < Aexp(-Bs). 

Proof. The representation of the Richardson model in terms of first passage perco- 
lation ensures the existence of A',B', M' > such that for each t > 0, 

(36) P(fjt £ B wt ) < A' exp{-B't). 

For more details, one can refer to Kesten [24] . 

We first control the process in integer times thanks to the following estimate: 

P(3k G N : m t B M >k+s/2) < G N : r)k+ s /(2M>) £ B M , k+s/2 ) 

+ 00 

k=0 

A' ( B's 

(37) " l-cxp(-B0 CXP r^ 

Let us now control the fluctuations between integer times. Let M > M': 
P{{3t > : m t B Mt +s} n {Vfc G N, Vk C B M , k+s/2 }) 

(38) < ]T F(3t G [fe, fc + 1] : ?7fe C B M , k+s/2 et r?* ^ B M t+a)- 

Then, denoting by C > a constant such that \B t \ < C'(l + t) d and by A, B the 
constants appearing in Lemma 131} 

P(3t G [fc, fc + 1] : 77 fe C B M , k+s/2 et ^ ^ B A/t+s ) 

< P(r/fc C B M , k+s/2 et ?7fc + i ^ B M k+ s ) 

(39) < |B M , fc+s/2 |P( m ^ S fe(M _ M , )+s/2 ) 

< C"(l + M'jfc + s/2) d Aexp(-B(k(M - M') + s/2)) 

< AC {I + s/2) d exp(-Bs/2){l + M'k) d exp(-B(k(M - M')). 



Inequality ([39]) comes from the Markov property and from the subadditivity of the 
contact process. Since the series ((1 + M'k) d exp(—B(k(Al — M')))k>\ converges, 
the desired result follows from 133 and 1381). □ □ 



6.2. A restart procedure. We will use here a so-called restart argument, which 
can be summed up as follows. We couple the system that we want to study (the 
strong system) with a system that it stochastically dominates (the weak system), 
and that is best understood. Then, we can transport some of the properties of the 
known system to the one we study: we let the processes simultaneously evolve and, 
each time the weaker dies and the stronger remains alive, we restart a copy of the 
weakest, coupled with the strongest again. Thus, either both processes die before 
we found any weak process surviving; in this case the control of large finite lifetimes 
for the weak can be transposed to the strongest one, or the strongest indefinitely 
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survives and is finally coupled with a weak surviving one. In that case, a bound for 
the time that is necessary to find a successful restart permits to transfer properties 
of the weak surviving process to the strong one. 

This technique is already old; that can be found for example in Durrett |14j . 
section 12, in a very pure form. It is also used by Durrett and Griffeath |16j . 
in order to transfer some controls for the one-dimensional contact process to the 
contact process in a larger dimension. We will use it here by coupling the contact 
process in inhomogeneous environment A <S A with the contact process with a 
constant birth rate A m i n . Here, the assumption A m i n > A c (Z d ) matters. 

To this end, we will couple collections of Poisson point processes. Fix A £ A. 
We can build a probability measure Pa on SI x fl under which 

• The first coordinate w is a collection ((w e ) egE d, {^ z ) z & d ) of Poisson point 
processes, with respective intensities (X e ) e(£E d for the bond-indexed pro- 
cesses, and intensity 1 for the site- indexed processes. 

• The second coordinate r\ is a collection {(j] e ) e ^d, (r/ z ) z( zz d ) of Poisson point 
processes, with intensity A m i n for the bond-indexed processes, and intensity 
1 for the site-indexed processes. 

• Site-indexed Poisson point processed (death times) coincide: for every z € 

z d , n z =lo z . 

• Bond-indexed Poisson point processed (birth-times candidates) are coupled: 
for each e £ E d , the support of r\ e is included in the support of u> e . 

We denote by £ A — £ A (ui,ri) the contact process in environment A starting from 
A and built from the Poisson process collection ui, and £ B = £ B (ui,r]) the contact 
process in environment A,„i n starting from B and built from the Poisson process 
collection r\. If B C A, then Pa almost surely, C,t C £ A holds for each t > 0. We 
can note that the process (£ , £ s ) is a Markov process. 
We introduce the lifetimes of both processes: 

t = inf{t > : Ct = I and. for x G z<i > T 'x = inf > : Cf = 0}- 

Note that the law of t' x under Pa is the law of t x under PA min ; it does actually not 
depend on the process starting point, because the model with constant birth rate 
is translation invariant. 

We recursively define a sequence of stopping times (iik)k>o and a sequence of 
points {zk)k>o, letting uq = 0, zq = 0, and for each k > 0: 

• if u k < +00 and £ Ufc 7^ 0, then u k+1 = r' Zk o Uk ; 

• if Uk = +00 or if £ Ufc = 0, then Uk+i = +00; 

• if Ufe+i < +00 and t; Uk+1 7^ 0, then z^+i is the smallest point of £ Uk+1 for 
the lexicographic order; 

• if itfc+i = +00 or if £, Uk+1 = 0, then Zk+i = +00. 

In other words, until Uk < +00 and £ Uk 7^ 0, we take in t; Uk the smallest point z k 
for the lexicographic order, and look at the lifetime of the weakest process, namely 
C, starting from z k at time Uk- The restart procedure can stop in two ways: either 
we find k such that < +00 and f Ufc = 0, which implies that that the strongest 
process (which contains the weak) precisely dies at time u k ; or we find k such that 
Uk < +00, £ Uk 7^ 0, and Uk+i = +00. In this case, we have found a point z k such 
that the weak process which starts from Zk at time Uk survives; particularly this 
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implies that the strongest also survives. We then define 

K = inf{?i > : it„+i = +00}. 

The name of the K variable is chosen by analogy with Section[3l The current section 
being independent from the rest of the article, confusion should not be possible. It 
comes from the preceding discussion that 

(40) (t = +00 ^=> ^ 0) and if r < +00, then u K = t. 

We regroup in the next lemma some estimates on the restart procedure that are 
necessary to prove Proposition [5j Recall that p is introduced in Equation ([7]) . 

Lemma 33. We work in the preceding frame. Then, 

• VA G A Vn G N V X {K > n) < (1 - p) n . 

• V-B G B{V) P A (r = +00, C ZK o 6 UK £B)=P A (r = +oo)P Amln (£° G B). 

• There exist a,/3 > smc/i that for every A G A, EA(exp(cm^ )) < /3. 

Proof. By the strong Markov property, we have 

P A (A">n+l) = P A K +1 <+oo) 

= PaK < +00, Cu„ ^ 0, <, o U J < +00) 

< V x (u n <+cx>)(l-p)=F x (K>n)(l-p). 

Thus, K has a subexponential tail, which proves the first point. Particularly, K is 
almost surely finite. 

Using lj40j) and the strong Markov property, we have also 

i\(r = +00, C K o 9 UK G £>) = P A (£„ K ^ 0, C" K o 0„ K G B) 

+OO 

-j-oc 

= EE < +(X3 ' C ^ ' ^ = z ' r L- ^ C ZK o UK G B) 

k=0 zGZ d 
+00 

= EE < C + ' Zfc = *) p v»> = e° g s) 

k=0 z<£Z d 

+00 

= PA mi „(r = +00, e G S)EPaK < +00, C- ^ )- 

fc=0 

Taking for i? the whole set of trajectories, we can identify: 

+00 

P(r = +00) = Pa(t = +c») = PA ml „(r = +c») ^P A («jk < +00,^ ^ 0), 

which gives us the second point. 

Since A m j n > A c (Z d ), The results by Durrett and Griffeath [16] for large A, 
extended to the whole supercritical regime by Bezuidenhout and Grimmett [4], 
ensures the existence of A, B > such that 

Vt>0 F Xtaln (t < t < +00) < Aexp(-Bt), 

which gives the existence of exponential moments for i~l{ r<+00 }. Since PA min (T = 
+00) > 0, we can chose (e.g. by dominated convergence) some a > such that 
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E A mi „(cxp(ar)l {r<+oo} ) =r < 1. 
For k > 0, we note 



Sk = cxp ^ r^. o 6» Ui j l{ Ufc <+oo}- 

We note that Sk is T Uh -measurable. Let k > 0. We have 

exp(au_ R -)l{ A - =fc} < Sk- 

Thus, applying the strong Markov property at time u^i < +oo, we get, for k > 1 

E A [exp(auK)l {K =fc}] < ^(5*^) = E A (S'fe_i)E Amta (exp(Q;r)l{ T<+00 }) 
< rE x (S k -x). 

Since r < 1, it comes that E A [exp(aitif)] < < +oo. □ 

6.3. Proof of Proposition [H Estimates ([8]) and (0 follow from a simple stochas- 
tic comparizon: 

Proof of ([7p. It suffices to note that for every environment AG A and each z G Z d , 
we have 



P A (r z = +oo) > Px mSa (r 2 = +oo) = P Amm (r° = +oo) > 0. 



□ 



Proof of (Hp. We use the stochastic domination of the contact process in environ- 
ment A by the Richardson model with parameter A max - For this model, (f36|) ensures 
a growth which is at least linear. □ 

Then, it remains to prove ((9j), (fTOj) and (fTTj) with a restart procedure. 

Proof of (G|). Let a, (3 > as given in the third point of Lemma [33l Recall that 
uk = t on {t < +00} . For each A G A and each t > 0, we have 

P A (t<r<+oo) = P A (e Qt <e QT ,r<+oo)=P A (e Qt <e QUK ,r<+oo) 
< P A (e Qt < e aUK ) < e~ at E x e aUK < /3e~ Qt , 

which concludes the proof. □ 

Proof of MQ\) . Since A m ; n > A c (Z d ), Durrett and Griffeath's results [16j for largeA, 
extended to the whole supercritical regime by Bezuidenhout and Grimmett [4], 
ensure the existence of constants A,B,c > such that, for each y G Z d , for each 
t > 0, 

(41) PA min (t(y) > M + 1 ) < Aexp(-Bt). 



Besides, the domination by the Richardson model with parameter A ma x and Lemmal32l 
ensure the existence of A, B, M > such that for every A G A, for each s > 0, 

(42) F x {3t > 0, &° t B M t+s) < AcM-Bs). 



ASYMPTOTIC SHAPE FOR THE CONTACT PROCESS IN RANDOM ENVIRONMENT 31 



By decreasing c or increasing M if necessary, we can also assume that -fj <1. Now 

Pa (t(y) > M + t , r 

tc\ - ( tc 

+P> in,- < 



+Pa (t = +c^, u K < C K C B tc/3 , t(y) > M + i 

The existence of exponential moments for m if given by Lemma [33] enables to bound 
the first term: there exist C, a > such that for each A € A, for each t > 0, 

~ / fc \ „ / act 



The second term is controlled with the help of (j42|) : 

«x < £ K t Btc/s) < Pa(3* > 0, # £ S Mt+f ) < Acxp (~Bj) • 

It remains to bound the last term. We note here 

t'{y) = inf{< > : y £ (?}■ 

Recall that if r = +oo, then £ UK ^ and zk is well-defined. Since t(y) is the 
hitting time of y and £° D Q for each t, we have, on {t = +oo}, 

t(y) < u K + t'(y - z K ) o T„ K o 9 UK . 

If ujr < ^ < |, then t(y) < | + i'fo - z^) o T ZA . o UK . If, moreover, ^ K C B te / 3j 
we have ||y|| > ||y — — which gives, with the second point in Lemma [33l 

Pa (t = +oo, UK < JjjL, e° K c S ct/3, t(y) >^+t 

< P A f t = +oo, t'(y - z K ) o T ZK o 6 UK > ^ ~* K ^ + | 

< P A (r = +oo) sup P Amin [ t(y - z) > t- — + | ] < Acxp(-B</2), 

zG z d V c 2J 

where the last inequality follows from (|4T]l . This concludes the proof. □ 

Proof of ill)) . Let s > 0, and denote by n the integer part of s. Let 7 > be a 
fixed number, whose precise value will be specified later. 



K' s ) = F(3t > s : i K t ) 
< H B ik <£K k ) + Y^ P{B lk C K k , 3t £ [fc, fc + 1) such that g K t ). 

k—n k—n 

Let us first bound the second sum. Fix fc > n. Assume that B lk C K k and 
consider t £ [fc, fc + 1) such that ^ A" t . Then, there exists x £ Z d such that 
£ £f\£t . Since £ £f and t > k, there exists y £ Z d such that y £ and 
£ £ t y _ fe o0 fc . If y £ B 7fe C AT fe , then = d (y) = 1, which implies that y £ ($.. 
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Now, since £ £f_j,°#fc, we obtain £ which contradicts the assumption ^ 
Thus, we necessary have y $ B^k, so: 

P\(B jk c K k , 3t e[k,k + 1) tel que £ iT t ) 
< -Pa(0g U ff\ B -rA 

p V se [°^] / 

P V "Si ' 1 ! / P 

Since the Richardson model with parameter A max stochastically dominates the con- 
tact process in environment A, we control the last term thanks to Lemma l3ll 

To control the first sum, it is sufficient to prove that there exist positive constants 
A, B, 7 - and this will fix the precise value of 7 - such that for each Ad A and each 
t > 

(43) Vx(B 7t <t K t , t° = +00) < Acxp(-Bt). 

The number of integer points in a ball being polynomial with respect to the 
radius, it is sufficient to prove that there exist some constants A, B, c! > such 
that for each t > 0, for each x £ Z d , 

(44) ||x|| < c't P A (£ t ° ^0,x£ < Aexp(-Bt). 

To prove i[44|) . we will use the following result, that has been obtained by Dur- 
rett [17] as a consequence of the Bezuidenhout-Grimmett construction [4]: if £ 
and £ x are two independent contact processes with parameter A > A c (Z d ), respec- 
tively starting from and from x, then there exist positive constants A,B,a such 
that for each t > and each x £ Z d , 

(45) \\x\\ <at=> Pfe° n £ = 0, |f ^ 0,£° ^ 0) < Acxp(-Bt). 

Let a and M be the constants respectively given by Equations 1(45]) and ((8|). We 
put c' = a/2 and choose e > such that c' + 2eM < a. 
Let a G B° at/4 and 6 G B» t/4 . We set 

a a ,s = C 0*t/2 and /? M = {y£Z d : b £ C v s o t(1 _ e/2) _ s }. 

Then, (a aiS ) < s < t / 2 (i_ e ) and (Ai,«)o<«<i/2(i-e) are independent contact processes 
with constant birth rate A m ; n , respectively starting from a and from b. The process 
(j3a,s)o<s<t/2(i-e) ls a contact process, but for which the time axis has been reverted. 
In the same way, we set 

g={y£Z d : xegoO t -.}. 
Note that (Cs)o<s<t/2 has the same law as (£?)o<s<t/2- Note that: 

• Assuming that a £ £° /2 , a a ,(i_ e ) 4 /2 n 0b,(i-e)t/2 ^ and & G /2 , then 
x£$. d 

• If x £ £f , then is non-empty. 

• If £° is non-empty, then £°, 2 is non-empty. 
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Thus, letting 

E° = {^^0}\{3«eB: t/4 nC o t/2 :« n ,( M(/2 ^0} 

and W = {if /2 ± 0} \ {3b e B x at/i fl |* /2 : ft i(1 - 8 )t/a ^ 0} , 

we get 

(46) p A (c t ° /0,xe #V) < pa(c? /2 + 0, i? /2 ^ 0, e t/2 n i? /2 = 0) 

< P A (£°) +Fa(^ x ) + S, 
where 5= ^ Pa (a^ d-t)- ^ 0, ^b. ' 1 -'" ^ ' a a , "7" n ' g b. (1 7 )t = )' 

For every couple (a, 6) that appears in S, we have \\a — b\\ < \\a\\ + \\b — x\\ + \\x\\ < 
at/4 + at/4 + at/2 = at, which allows to use lj45|). and gives the existence of 
constants A,B,C > such that 

S < C"(l + at/4) 2d Aexp{-B(l - e)t/2). 

By another time reversal, we see that Pa^ 1 ) = F X .\(E°); then it suffices to control 
¥\(E°) uniformly in A. Let 

Ei = {it/2 + 0}\{3a eZ d : a G £° /2 , a a , ( i_ e )t/ 2 ^ 0}- 

We have P A (£°) < Pa(#i) + i\(£° t/2 ^ B °t/ 4 )- % the choice we made for £ and 
Inequality (HJ), we have 

VAeA Vt>0 F x (g t/2 <£B(0, at/4)) < Aexp(-Bst/2). 

Thanks to the restart Lemma l33l we can see that 

f>x(u K > et/2) < (3exp(-ast/2). 

Suppose then that uk < et/2 and £®, 2 ^ 0: zk is thus well defined and we have 
t' Zk o 9 Uk = +oo. Then, there exists an infinite infection branch in the coupled 
process in environment A m i n starting from £® K . This branch contains at least one 
point a <E £[i_ £ ) t / 2 - By construction a £ £(i_ E w 2 an d a a (i_ E ) t / 2 7^ 0, which 
completes the proof of lf43|) . □ 

Remark On our way, we proved that for each A 6 A, 

lim P A (£° ^ 0, Q * 0, $ n |f = 0) = 0, 

t — >-t-oo 

which is the essential ingredient in the proof of the complete convergence Theorem[4j 
One can refer to the article by Durrett [17] for the details in the case of the classical 
contact process. 

Appendix: proof of almost subadditive ergodic Theorems 1231 and l24l 

Proof of Theorem\23[ Let a p = C p ^ a and u n — E[/ n ]: for every n,p <E N, we have 
E[r+ p ] < (n(r+ P ) a ]) 1/a < Cl /a = a p , hence 

u n+p < u n + u p + E[g p ] + E[r n>J) ] < u n + u p + E[g p ] + a p . 
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The general term of a convergent series tends to 0, so C p = o(p a ), or a p = o(p). 
Since a " +Eg " tends to 0, the convergence of u n /n is classical (see Derriennic [TT] 
for instance). The limit [i is finite because u n > cn holds for each n. 

We are going to show that / = lim ^ stochastically dominates a random 

n — >+oo 

variable whose mean value is not less than fx. 

For every random variable X, let us denote by C(X) its law under P. We denote 
by /C the set of probability measures on whose marginals m satisfy: 

Vt >0 m(]£,+oo[) <P(/i+5i > t/2) + d{2/t) a . 

Define, for k > 1, 

Afc = /fe+i - A-, 

and denote by A the process A = (Afc)fc>i. For k £ N, subadditivity ensures that 
Afc < (fi + gi) o 9 k + r k ,i, hence, for each t > 0, 

P(A fc >t) < P((/i+ 5 i)ofl fc >i/2)+P(r+ 1 >t/2) 
< P(/i+ 5 i >t/2) + Ci(2/*) Q - 

This ensures that A G /C. 

We denote by s the shift operator: s((ufe)fe>o) = (u/c)fe>i, and consider the 
sequence of probability measures on K N : 



{L n )n>X = (i^£(s^oA) I 

V^ 1 I n>l 



Since K, is convex and invariant by s, the sequence (£ n ) n >i is A^-valued. Let 
n, k > 1. 

„ n 
TT k (x)dL n (x) = -^E(7Tfc(^oA)) 

71 3=1 

n 

- £E(/ fc+j+1 - A+,) = -(E[/ n+fc+1 ] - E[/ fc+1 ]). 



n 

3=1 



Let M fc = sup I|E[/„ +fe+ i] - E[/ fe+ i]|. The convergence of u„/n implies that M k 

n>l 

is finite. Similarly, the subadditivity gives 

r l - 

n+(x)dL n (x) = -^E(7r+(^oA)) 
i=i 
1 " 



n 

3=1 

Thus, we have 

|7r fc (x)| dL n {x) < I 2tt+(x) dL n (x) + \ J n k (x) dL n (x)\ 

< M k + 2E[f+] + 2E[g+] + 2a x . 

Let K! be the family of laws m on M N such that for each k, J \n k \ dm < 2M k + 
+^[gt] +Oi . KJ is compact for the topology of the convergence in law and the 
sequence {L n ) n >\ is /C'-valued. So, let 7 be a limit point of (L n ) n >i and {nk)k>i 
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a sequence of indexes such that L nk =>■ 7. By construction, 7 is invariant under 
the shift s. 

Now, the sequence of the laws of the first coordinate wi(x) under (L nk )k>o 
weakly converges to the law of the first coordinate under 7. Also, by definition 
of /C, the positive parts of these elements form an uniformly integrable collection, 
so fir* = lim J tt^ dL nk . However the Fatou Lemma tells us that / 7rf dj < 

lim J ir^ dL Uk , hence finally 

k — >+oo 



7Ti c?7 > lim / 7Ti dL nk = fl. 

k— >+oo J 

Let Y = (Yfe)fe>i be a process whose law is 7. Since 7 is invariant under the 
shift s, the Birkhoff Theorem tells us that the sequence (i X)fc=i ^fe)n>i a - s - con- 
verges to a random variable Y x , which satisfies then E(Yoo) = J tti > fx. 

It remains to see that the law of Yoo is stochastically dominated by the law of 

/ = lim We will show that for each a G M, P(F OCl > a) < P(/ > a). By 

n — >+oo 

left-continuity, it is sufficient to prove the inequality in a dense subset of R. Thus, 
we can assume that a is not an atom for the law of /. 

{Yoo > a \ = < hm > a > = (J < mi > a 

I n J fc>i l_ «>fe n 

Hence 

P(Foo > a) 

= EE IV f mf 7ri + -" + 7r " > a 

lim" inf P y ( inf 7ri + 7 + 771 > a 

k^+00 n>k \ k<i<n I 

1 nK ( 4 4 

< Hm inf lim Vp inf — ZllLZll s j a > 

fe^+oo n>fc if ^+00 riif V k<i<n i 

3=1 



Let e > 0. We have, for fixed k, n,j: 



. , TTl H h 7Tj ■ 

mi ; o s J o A > a 

k<i<n i 

inf >fl 

k<i<n i 

inf CA±M^±l±l2±M > 

k<i<n I 

inf (/l+9l)0 ^>a- £ Vp( sup!>LM>; 

k<i<n i I \ i>k % . 
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On one hand, we have 

r j+i,i 



sup 

i>k i 



>e < 



E 

i>k 



r+ 



>e a < e- E ^E[(r; +M )«] 



e - a y 9l 



We can note that this term does not depend on j nor on n. On the other hand, 



• r ih + 9i) ° Oj+i . 
mi — > a 



k<i<r 



■ t fi+9i . 

mi ■ — > a — e 



k<i<n i 



which does not depend on j. Then, for each e > 0, we have for every n,k, with 
n > k: 

lim —Yp( inf 7ri + ' ' ' + ^ o J o A > a 

m.i» ^ — ' V k<i<n i 



A'^+oo Tiff 



£ 



. , ft+9i ^ 

mi ■ > a — e 



k<i<n I 



next 



1 nK / 
inf lim — Vp inf 



< 



n>k K— >+oo 71if 



7Tl 



— o s- 7 o A > a 



k<i<n 



V _L + inf P inf 

i>i i Q n>fc V k<i<n 



fi + 9i 



> a — e . 



Finally, 



P(F 00 > a) < lim inf 



• f /» + ft 

ml ■ > a — £ 



fe — >+oo n>k V k<i<n % 

f >a-£ 



ImT e" a E 



fc->+oo V i ~^- k 1 



< 



lim 



fi + 9 1 



i — >+oo I 



>a-£|=P( lim — >a-£ 



considering that gi/i almost surely converges to 0. Letting e tend to zero, we obtain 
P(F 00 > a) < P [ hm — > a] = F{f > a). 
It remains to see that / is invariant under the 6 n 's. Fix n > 1. We have 



+ °° /r+ 
/ ' n,p 

p— 1 v 



n,p 

p 



Particularly, - J ^ £ almost surely converges to when p tends to infinity. Since 
fn+p < fn + fp ° 8 n + g p o 9 n + r+p, dividing by n + p and letting p tend to +oo, it 
comes that 

£<f_°8 n a.s. 

Since P is invariant under 9 n , we classically conclude that / is invariant under 9 n . 

□ 
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Remarks. In the present article, we made no use of the possibility to take a non- 
zero g p . In the case where the (g p ) are not zero, but the r n ,p'& are, we obtain a result 
which sounds much like Theorem 3 in Schiirger [34]. Like Schiirger [33], we use the 
idea of a coupling with a stationarized process. This idea is due to Durrett [13] and 
has been popularized by Liggett [28]. However, there is here a refinement, because 
we directly establish a stochastic comparison with the random variable Y, whereas 
previous papers establish a stochastic comparisons with the whole process {Y n ) n >i, 
that admits Y as its infimum limit. 

In most almost subadditive ergodic theorems, almost sure convergence requires 
strong conditions on the lack of subadditivity (stationarity for instance) . We obtain 
here an almost sure behavior by only considering a condition on the moments (of 
order greater than 1) of the lack of subadditivity. Besides, we know that bounding 
the first moment of the lack of subadditivity is not sufficient to get an almost sure 
behavior (see the remark by Derriennic [11] and the counter-example by Derriennic 
and Hachem |12j). 

Proof of Theorem\2^ It remains to prove that E ( Em —J < /i. 

\ n — >+oo Tl J 

We fix k > 1. By subadditivity, we have for each n > and every < r < k — 1: 

fnk+r < fnk + (fr+gr)°Onk+r^ kr 



< 



Since P is invariant under 
convergence: 




fnk — fk ° 
^ i=0 



;)<] +(fr + 9r)°e nk + r+ 



the Birkhoff Theorem gives the L 1 and almost-sure 



lim 



iv 



fk o 



where Xk is the cr-algebra of the ^-invariant events. Let us now control the resid- 
ual terms. Since the finite collection {f r + g r )o< r <k-i is equi-integrable and P is 

invariant under 8k, the collection ( sup (f + q ) o 0?) n>1 is equi-integrable, 

0<r<fc-l 

which ensures the almost sure and L 1 convergence: 



1 



lim — sup (f r 

n— +oo n <r<k-l 



+ g r )o(0 k ) n =0. 



We have ^ E 



r+ 

^ nk,r ^ a 



< y ^ < 



-oo, which implies, as previously, that 



nk,r 



hence 



/n almost surely converges to 0. Finally, 



Vr e {0, 



lim — 



< 



nf k 



.fa 



lim - '" k+r 
; ^ +00 nk + r 



< 



E[f k \l k 



. We complete the proof by letting k tend to +oo. 

□ 



Remark. When there is no lack of subadditivity, the assumptions of Theorem l24l 
obviously hold; we obtain thus a subadditive ergodic theorem which sounds very 
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much like Liggett's |28]. However, these theorems are not strictly comparable, in 
the following sense that no one implies the other one. 

Indeed, extending a remark made by Kingman in his Saint-Flour's course [26] 
(page 178), we can note that the assumption of Kingman's original article - the 
stationarity of the doubly indexed process (X Si t) s >o,t>o _ can be weakened in two 
different ways: 

• Either assuming that for each k, the process (^( r -i)fc,rfc)r>i is stationary 
- this assumption will be used by Liggett [28] . 

• Or assuming that the law of X n , n + P does not depend on p. That assump- 
tion, suggested by Hammersley and Welsh, is the one that we use here, also 
used by Schiirger in |34j . 

Note however that the special assumption of stationarity is used in Liggett's proof [28] 
only in the so-called easy part, i.e. the bound for the supremum limit. 

Kingman thought that the first set of assumptions surpassed the second one, in 
view of possible applications. More than 30 years later, the progresses of subadditive 
ergodic theorems, particularly about bounding the infimum limit, lead to moderate 
this affirmation. 
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